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A  macroscopic  description  of  the  linear  Weibel  instability,  based  on  a  closed  set  of  linear  moment 
equations,  is  presented.  The  moment  equations  are  derived  from  the  linearized  Vlasov  equation  by 
taking  the  appropriate  velocity  moments  of  it  and  the  closure  is  achieved  by  means  of  an 
assumption,  which  is  justified  when  the  temperature  anisotropy  is  strong.  The  macroscopic 
description  is  manifestly  more  informative  of  the  physical  mechanism  of  the  instability  than  the 
kinetic  description.  It  is  hoped  that  the  researchers  will  find  such  a  description  analytically  more 
convenient  to  use  in  solving  plasma  physics  problems  where  Weibel  instability  due  to  strong 
temperature  anisotropy  plays  a  role.  [DOI:  10.1063/1.1521716] 


Particle  distributions  in  velocity  space  with  an  enhanced 
temperature  along  some  direction  are  quite  common  in  both 
laboratory  and  space  plasmas  where  the  degree  of  colhsion- 
ality  is  relatively  low.  There  are  a  number  of  ways  in  which 
such  a  temperature  anisotropy  can  develop  in  magnetically 
confined  as  well  as  in  magnetic  field-free  plasma.^  Coulomb 
collisions  eventually  make  the  plasma  distribution  isotropic. 
However,  collective  processes  (instabihties)  due  to  the  tem¬ 
perature  anisotropy  can  be  more  effective  than  the  binary 
collisions  in  moving  the  plasma  toward  an  isotropic  state. 
These  instabilities  can  be  both  electrostatic  and  electromag¬ 
netic  in  nature.^  Here  we  shall  deal  with  a  particular  electro¬ 
magnetic  instability  due  to  temperature  anisotropy,  known  as 
the  Weibel^  instabiUty,  which  is  excited  in  collisionless 
plasma,  even  in  the  absence  of  an  external  magnetic  field. 
The  unstable  waves  are  transverse  electromagnetic  waves, 
involve  only  the  electron  population,  and  do  not  produce 
density  perturbation.  The  Weibel  instability  is  considered  to 
be  one  of  the  mechanisms  for  magnetic  field  generation^’^  in 
laser  produced  plasma,  with  special  significance  in  laser  fu¬ 
sion  experiments.  It  may  also  play  an  important  role  in  the 
magnetic  reconnection  process  in  the  Earth’s  magnetotail  by 
participating  in  the  dynamics  of  the  electrons  in  the  so-called 
electron  diffusion  region.^ 

Weibel’s  work  stimulated  a  series  of  further 
investigations^"^"^  of  the  transverse  electromagnetic  instabil¬ 
ity  in  unmagnetized  plasma.  These  papers  dealt  with  the  lin¬ 
ear,  quasilinear  and  fiilly  nonlinear  theories  as  well  as  the 
computer  simulation  experiments  of  the  instability.  At  the 
same  time,  several  other  authors^^"^^  investigated  the  elec¬ 
tromagnetic  instabilities  in  magnetized  plasma  for  a  wide 
variety  of  anisotropic  velocity  distributions  and  for  different 
orientations  of  the  propagation  vector.  Most  recently,  Cali- 
fano  have  carried  out  further  investigations  of 

Weibel-type  instability,  where  the  role  of  temperature  anisot¬ 
ropy  is  taken  by  two  counterstreaming  electron  populations. 
All  of  the  previous  analyses,  including  that  of  Weibel,  have 
been  based  on  the  Vlasov—Maxwell  formalism. 

“^Electronic  mail:  bainandas.basu@hanscom.af.rml 


In  the  Vlasov-Maxwell  formalism,  the  linear  dispersion 
relation  for  Weibel  instability  in  unmagnetized  plasma  is  de¬ 
rived  by  solving  the  linearized,  nonrelativistic  Vlasov  equa¬ 
tion, 

(|.V.  v)7(r.,.,)- ^(e+ iTXB) . 

and  the  two  Maxwell’s  equations, 

-  1 

VXE=--— ,  (2) 

c  dt 

_  4’7re  f  _  1 

VXB= - —j  dv\f(r,v,t)+-—,  (3) 

self-consistently.  Here  /o(v)  is  the  unperturbed  electron  ve¬ 
locity  distribution,  7(r,v,r)  is  the  perturbation  of  the  distri¬ 
bution,  and  E  and  B  are  the  perturbed  electric  and  magnetic 
fields,  respectively.  Considering 

no  I  mV'^  \  m  ,  ,  m  , 

(4) 

and  space-time  dependence  of  the  perturbations  of  the  form 
~exp[/(fe““wO]»  the  linear  dispersion  relation  of  the  trans¬ 
verse  waves  (k»E=0)  is  then  obtained  as 

l  +  ^0.  (5) 

Here  =  is  the  electron  plasma  frequency, 

y7’^  =  (2r||/m)^^  is  the  electron  thermal  speed  associated 
with  r,, ,  and  W(^)  =  - 1  -  ^Z(^),  where  Z(|)  is  the  plasma 
dispersion  function.  WeibeP  demonstrated  instability  by  con¬ 
sidering  the  \  (jD/{kVT^^\>l  limit  of  Eq.  (5).  A  more  system¬ 
atic  analysis  of  Eq.  (5)  can  be  found  in  Ref.  14.  The  results 
are  summarized  below. 

For  small  deviations  from  isotropy,  (r^““r||)/r||<^l,  in¬ 
stability  is  found  for  a  low  frequency  mode  such  that 
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|a)/(A:V7’^)|<^l.  With  the  substitution  of  1 

for  and  the  use  of  a><^c/:,  Eq.  (5)  becomes 


,  ^pe  T'lI  <0 

'"n 


where  k>0  has  been  assumed.  The  solution  for  <a  is  then 
given  by'^ 

It  shows  that  purely  growing  waves  are  excited  when  Tj_ 
>r||  and  for  the  values  of  k  in  the  range  0<k<k„ ,  where 
kl={io)jc^){TJT^-\). 

For  strong  anisotropy,  (rx-r||)/r||>l,  the  approxima¬ 
tion  |a»/(fcV7-^)|<^l  with  ft)  given  by  Eq.  (7)  tends  to  break 
down,  except  when  k<^k„.  For  k<k„,  the  instability  is 
found  for  frequency  ft)  such  that  |ft)/(kV7-^)|^l.  In  this  re¬ 
gime,  and  so  Eq.  (5)  becomes 


ft)2-c2k2- 


For  <a<ck,  this  again  yields  a  purely  growing  solution  given 

by' 


0}=ik\  — 
\  m 


It  may  be  verified  that  the  restriction  |<o/(itV7-  )|>  1,  with  o) 
given  by  Eq.  (9),  is  satisfied  only  when  . 

In  this  Brief  Communication,  we  present  a  macroscopic 
description  of  the  linear  Weibel  instability  for  the  strong  an¬ 
isotropy  (7j^>r|j)  case.  The  description  is  based  upon  a 
closed  set  of  linear  moment  equations  that  are  derived  from 
the  linearized,  nonrelativistic  Vlasov  equation  by  taking  the 
appropriate  velocity  moments  of  it.  The  derivation  requires 
integration  by  parts  and  use  of  the  synunetry  properties  of /o 
in  velocity  space.  The  closure  is  achieved  by  means  of  an 
assumption,  which  is  justified  when  \(o/(kVr^^\>l,  i.e., 
when  the  temperature  anisotropy  is  strong.  A  similar  macro¬ 
scopic  description  for  the  weak  anisotropy  case  is  not  pos¬ 
sible  since  the  instability  relies  on  a  wave-particle  resonance 
process. 

We  consider  transverse  waves  propagating  in  the 
z-direction,  so  that  By  taking  the  successive  ve¬ 

locity  moments  of  Eq.  (1),  we  find  the  following  intercon¬ 
nected  chain  of  moment  equations: 

^  I  d  ^  e  ^ 

[(BXPo)„+(BXPo)  J.  (11) 

^  ^  f  e  _ 

‘/vv^V(r,v,r)  =  --Pozz£z.  (12) 

and  so  on.  These  are  the  linearized,  one-dimensional  versions 
of  the  moment  equations^'*  that  can  be  obtained  from  the 


exact  Vlasov  equation.  Here  is  the  x-component  of  the 
perturbed  electron  mean  velocity  u,  is  the  xz-element  of 
the  perturbed  stress  tensor  P,  is  the  xzz-element  of  the 
perturbed  heat  flow  tensor  Q,  and  their  standard  definitions 


Uxir,t)=—  f  dwJ{T,v,t), 
fto  J 


^xz(’"’0  =  m  j  dw^vj(r,v,0, 

Qxzz(*‘’0  =  mj  dw,v^7(r,v,r). 


The  unperturbed  stress  tensor  Pq  ,  corresponding  to  the  /q 
given  by  Eq.  (4),  is  diagonal  and  its  elements  are  Pqjcx 
~^Oyy—^o^i  Other  elements  being 

zero.  The  notations  on  the  right  hand  side  of  Eq.  (11)  mean 
xz-  and  ZA:-elements  of  the  tensors  resulting  from  the  cross 
products.  The  moment  equations  are  the  exact  consequences 
of  taking  the  appropriate  velocity  moments  of  Eq.  (1),  using 
the  symmetry  properties  of/o  in  velocity  space.  However, 
they  do  not  form  a  closed  set,  since  each  equation  contains  a 
term  of  the  higher-order  moment.  We  truncate  the  chain 
of  equations  by  assuming  m(<9/<?z)/t/w^v^7(r,v,r) 
so  that  Eq.  (12)  reduces  to 

~  e 

^^jczz  ~  ^  ^Ozz^x  •  (16) 

Simple  order-of-magnitude  estimates  indicate  that  the  valid¬ 
ity  of  the  assumption  requires  \(of{kVp^)\>l  for  a  wave 
propagating  along  the  z-direction  with  frequency  (o  and  wave 
number  k.  We  shall  examine  this  assumption  in  more  details 
later  in  the  paper. 

Equations  (10),  (11),  and  (16)  together  with  the  two 
Maxwell’s  equations,  c{VXE)=  -  dB/dt  and  c(VXB) 
=  —  4'n‘enQn-\-  SE/ dty  provide  a  closed  macroscopic  descrip¬ 
tion  of  the  linear  Weibel  instability  in  the  regime  correspond¬ 
ing  to  \o>/(kVY^)\>l,  To  verify  this,  we  substitute  solutions 

of  the  form  A(T,t)=Aik,o))exp[i(kz-(ot)]  into  Eqs.  (10), 
(11),  and  (16)  to  find 


2  = - p - p 

mn^d)  ^  mo)  ^ 


.  k  ^  ieriQ 
'’«=:;;2.zz+— (ri-rjB,, 


Qxzz=- 


From  Eqs.  (18)  and  (19)  we  have 

.  icmo^T’ii  ..  ieiiQ 

Pxz= - h^Ex+ - ^(7-||-7’^)B, , 

‘  mo)^  meat 

and  then  combining  Eq.  (20)  with  Eq.  (17)  we  find 


,  ,  *(7’l|-7’i)  J 

IH - T  Er - By,  . 

mo)  mc(o  ^ 
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Finally,  combining  Ej^—{a)lck)By  and  By— —  {Am! 
ck)enf^UxA‘{(j)lck)Ey.,  which  follow  from  the  Maxwell’s 
equations,  with  Eq.  (21)  we  obtain  the  dispersion  relation 

1  + - 7]  =0,  (22) 

which  is  exactly  the  dispersion  relation  that  is  obtained  from 
the  Vlasov-Maxwell  formalism  in  the  limit  \<ol{kVT^\>\ 
[see  Eq.  (8)].  We  recall  that  the  condition  \(ol{kyT^\>\ 
corresponds  to  the  strong  temperature  anisotropy  case. 

Let  us  now  return  to  the  assumption  by  which  the  term 
m{dl dz)i dyyy^^J{T,y,t)  in  Eq.  (12)  was  neglected  and  thus 
closure  was  achieved.  If  we  multiply  Eq.  (1)  by  ,  inte¬ 
grate  over  velocity  space  using  the  symmetry  properties  of /o 
and  neglect  the  higher-order  velocity  moment  term,  we  find 

d  f  -I*.  2>€TIqT\\  ^ 

-J  (23) 

So,  had  we  retained  the  term  m{dldz)idWj^v\f{Tyyyt)  in  the 
equation  for  then  instead  of  Eq.  (20)  we  would  have 
obtained 


Pxz=- 


mo)  mco)  ^ 


+ 


ieriQ 

mco) 


mo)^  j 


(r„-rj5,. 


(24) 


where  the  last  term  on  the  right  hand  side  of  Eq.  (24)  is  the 
additional  term.  It  is  clear  that  the  additional  term  may  be 
neglected  in  comparison  with  the  second  term  when 
\o)/{kVT^^\>l,  Similarly,  it  can  be  shown  that  the  contribu¬ 
tions  to  higher-order  velocity  moment  terms  are  of 

higher  order  in  e=|^y7’^^/ft>|<^l,  and  so  are  even  smaller. 
Hence,  the  closure  of  the  chain  of  moment  equations  by 
adopting  Eq.  (16)  is  indeed  justified  when  \a}/(kVj^^\>l, 

It  should  be  pointed  out  that  an  equivalent  macroscopic 
description  of  the  instability  is  obtained  if  the  equations  for 
Uy ,  Pyz  and  Qyzz  are  considered  instead.  These  equations  are 
readily  obtained  from  Eqs.  (10),  (11),  and  (16)  with  x  re¬ 
placed  by  y  everywhere.  In  other  words,  there  are  two  lin¬ 
early  independent  eigenstates,  represented  by  {Ux  ^Pxz^Qxzz} 
and  {uy.PyzjQyzz}^  respectively,  with  the  same  eigenvalue 
(Ok .  This  degeneracy,  which  is  a  consequence  of  the  symme¬ 
try  of  /o(v)  in  velocity  space,  is  removed  when  an  external 
magnetic  field  Bq  is  applied.  The  moment  equation  descrip¬ 
tion  presented  above  can  be  extended  to  the  magnetized 
plasma  case  in  a  straightforward  manner  by  noting  that  Eq. 
(1)  will  have  an  additional  term  represented  by  —{el me) 
X(vXBo)*(^7/^v),  and  can  be  used  to  study  the  effects  of 
the  external  magnetic  field  on  the  Weibel  instability.  For  ex¬ 
ample,  if  we  take  Bq  to  be  along  the  z-direction,  the  addi¬ 
tional  term  leads  to  coupling  between  the  states 

{ux^PxzMxzz}  {^y^Pyz^Qyzz}  cyclotron  motion 

of  electrons.  Consequently,  the  corresponding  equations  for 
Uy ,  Py^  and  Qy^^  are  needed.  Those  equations  can  be  derived 
in  a  similar  manner  using  the  same  closure  approximation.  It 


is  found  that,  in  magnetized  plasma,  circularly  polarized 
waves  are  the  proper  eigenmodes  with  different  dispersion 
relations  for  the  right  and  the  left  circular  polarizations.  Fur¬ 
thermore,  the  perturbed  electric  field  associated  with  the 
transverse  waves  propagating  only  along  the  z-direction  has 
no  electrostatic  part.  It  arises  solely  from  the  fluctuating 
magnetic  field. 

The  macroscopic  description,  based  on  the  moment 
equations,  is  manifestly  more  informative  of  the  physical 
mechanism  of  the  Weibel  instability  than  its  kinetic  descrip¬ 
tion.  The  proposed  moment  equations  and  the  subsequent 
analysis  leading  to  the  derivation  of  the  dispersion  relation 
show  that  the  stress  tensor  term  P^z  (or  Py^)  in  the  momen¬ 
tum  balance  equation  plays  a  crucial  role  in  the  excitation  of 
the  instability.  Without  this  term,  one  would  recover  the  fa¬ 
miliar  stable  electromagnetic  mode  in  plasma,  which  is  de¬ 
scribed  by  (o^  —  c^k^-\-(Opg .  The  analysis  further  shows  that, 
in  order  to  accurately  determine  P^z  (or  Py^),  the  heat  flow 
term  Q^zz  (or  Qyzz)  describing  perpendicular  (to  k)  flow  of 
the  pardlel  (to  k)  thermal  energy  due  to  the  transverse  elec¬ 
tric  field  Ex  (or  Ey)  must  be  retained.  The  physical  mecha¬ 
nism  represented  by  the  moment  equations  is  that  the  force 
experienced  by  the  electrons  in  the  fluctuating  field  and  the 
associated  heat  flux  cause  a  momentum  flux  (P),  which  af¬ 
fects  u  (and  hence  current  density  j)  in  such  a  way  as  to 
increase  the  field  fluctuation.  Earlier,  Fried^  offered  a  simi¬ 
lar  explanation  by  means  of  an  approximate  treatment,  and 
the  role  of  heat  flow  was  not  recognized.  The  moment  equa¬ 
tion  description  presented  here  seems  to  provide  a  more 
complete  and  accurate  picture  of  the  instability.  A  kinetic 
description  of  linear  Weibel  instability  can  become  quite 
cumbersome  analytically,  if  the  equilibrium  magnetic  field  is 
inhomogeneous  as,  for  example,  is  the  case  with  the  Earth’s 
magnetotail,  since  it  involves  the  calculation  of  unperturbed 
particle  orbits  and  the  integration  of  the  linearized  Vlasov 
equation  along  those  orbits.  The  same  may  be  true  for  other 
space  and  laboratory  plasmas.  So,  it  is  hoped  that  the  re¬ 
searchers  will  find  the  moment  equation  description  analyti¬ 
cally  more  convenient  to  use  in  solving  plasma  physics  prob¬ 
lems  where  Weibel  instability  plays  a  role.  But,  we  reiterate 
that  the  applicability  of  this  description  is  limited  to  the  hy¬ 
drodynamic  regime  [\(o/{kVf^^\>l],  which  in  the  context 
of  the  Weibel  instability  corresponds  to  strong  temperature 
anisotropy. 
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